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Abstract 

We present integrable lattice equations on a two dimensional square lattice with 
coupled vertex and bond variables. In some of the models the vertex dynamics is 
independent of the evolution of the bond variables, and one can write the equations as 
non-autonomous "Yang-Baxter maps". We also present a model in which the vertex 
and bond variables are fully coupled. Integrability is tested with algebraic entropy as 
well as multidimensional consistency. 



1 Introduction 

\ Along with the classification results for integrable lattice maps on quad-graphs with one 
component fields on vertices [T], various more general integrable lattice equations have been 
proposed recently. This includes the "Yang-Baxter maps" (solutions to the functional Yang- 

'. Baxter equations [2]) with variables on bonds [21 HI [5], and number of higher order as well 

', as multi-component cases [HI [3 [8] . 

We describe here integrable lattice models which have both vertex and bond variables^ 
These models were found in the analysis of a coarse graining process of known integrable 
models on a square lattice [10] . We first describe the basic procedure we used, and the way to 
test integrability via algebraic entropy (section |2]) . We then list various models obtained in 

*E-mail: jarmo.hietarinta@utu.fi 
^E-mail: claude.viallet@upnic.fr 
"'"permanent address 

^ The term "lattice models with bond and vertex variables" has already been used in the context of 
quantum spin models [5], but the setting is then quite different: the vertex and bond variables are independent 
variables and take on discrete values, the main data being the Boltzmann weight of the configuration, and 
one is the interested in the statistical mechanics of the system. Here, in contrast, the lattice equations define 
a 1+1 dimensional evolution of the variables. 
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this way, giving explicitly the evolution they define and their entropy. We also show that the 
evolution may be determined from a set of algebraic relations which have special rationality 
properties (section [3]). We finally examine the three dimensional consistency of these models, 
and show that some of them can be interpreted as solutions of non-autonomous functional 
Yang-Baxter equations (section Hj). 

2 A road map 

2.1 The setting 

The starting point is a regular planar square lattice. Its vertices are labeled by two integers 
{n,m). We introduce two kinds of variables: the vertex variables, denoted Wn,m, and the 
bond variables denoted Xn^m for the horizontal bonds, and Fn,m for the vertical bonds. The 
bond variables are indexed by their base vertex, which is to the left (respectively down) 
along the bond, as shown in Figure [U^a). As usual we use the shorthand notation where only 
shifts are indicated, e.g., X2 = Xn,m+i,wi2 = Wn+i,m+i, etc... 

The models are given by three relations between the eight variables associated to each 
square cell, in such a way that they define an evolution: from the three defining relations, 
it is possible to calculate X2, Yi, W12 (open circles in Figure [TJ^a)) from w, X, Y, wi, W2 (black 
disks). To fully define an evolution, we have to give a suitable set of initial data. Such data 
may be given on a diagonal staircase, after which the vertex and bond variables may be 
evaluated on the entire lattice, using the local relations on each cell. 

We will be interested in models where the evolution is rational, i.e., X2,Yi,wi2 are ra- 
tionally expressed in terms of w, X,Y,wi,W2. For the models we will consider, one may 
actually permute the role of the four corners, and in particular define a backward evolution 
calculating w,X,Y in terms of Yi, X2,wi,W2,wi2- This is a form of rational invertibility, 
similar to birationality for maps. For the models we describe, it is actually possible to define 
rational evolution in all four directions on the lattice. 

2.2 A construction procedure 

There is one very simple way to generate such integrable maps from known one-component 
lattice maps. 

Suppose we are given an integrable quad map. As shown in Figure HJ^b), consider 
{a;, xii, X22, 3^1122} as vertex variables of a coarser lattice (renaming them {w,wi,W2,wi2}, 
the rule being Wn,m = X2n,2m), and {x2, Xu2, xi, X122} as bond variables (and renaming them 
{Y,Yi,X,X2}, i.e., Xn,m = X2n+i,2m, Yn,m = X2n,2m+i)- The quad map determines Xu in 
terms of x, Xi,X2, and 2:122 in terms of X2, X22, X12, and so on, after which one may forget the 
central X12. One signature of this construction is that X2 does not depend on wi, and Yi 
does not depend on W2- 

The purpose of this paper is to present a number of models which cannot be obtained 
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Figure 1: (a): the basic lattice square, and (b): its precursor, a 2 x 2 sublattice. 



in the simple way described above, and to analyze their properties. They will in particular 
satisfy 

- 7^ 0, Vn, m. [1) 
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In order to construct such maps we start from quad maps having a 2 x 2 Lax pair [TT|[T2]. 
We use this Lax pair to write the zero curvature condition on the larger square of Figure 
[UJ^b). From the size of the Lax matrices, we see that this zero curvature condition can give 
at most three scalar conditions. We know that one rational solution can be obtained in the 
way described above, we call this solution "regular" . What is remarkable is that for certain 
quad maps, there exists another rational solution, which we called "exotic" |10] . 

In other words, we use the fact that the zero curvature condition written on a coarser 
lattice is ambiguous. In most cases there is only one solution giving rational evolution, it 
is the regular one. There are usually other solutions that are not rational, but we discard 
them. Only in some exceptional cases do we find rational solutions satisfying 



2.3 Testing integrability 

When an exotic solution exists, there is no guarantee a priori that it leads to an integrable 
evolution, and we have to check its integrability. 

The setting is appropriate to use the vanishing algebraic entropy criterion [HKTS]. Given 
initial data on a basic diagonal staircase, one evaluates the w's, X's and F' s away from the 
initial diagonal in terms of these data. The entropy can then be extracted from the sequences 
{dn} of the successive degrees respectively for the w's, X's and F's. The entropy is the limit 

?7 = lim - \og{dn). 

n— !>oo n 

This limit always exists. The characteristic of integrability is the vanishing of f]. It is 

^In some respects this is similar to using conservation laws to derive dynamics, which sometimes allows 
multiple solutions jl3j. 
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equivalent to polynomial growth of the sequences {dn}- Non- vanishing of the entropy, that 
is to say exponential growth of a sequence {dn}, means non-integrability. 

Remark : The simplest way to find the exact value of the entropy is by fitting, if possible, 
the generating function ({s) = J2n^nS^ sequence of degrees by a rational fraction. 

The entropy is then determined by the position of the poles of that fraction [HI [15] . 

2.4 Simplified defining equations 

We get the exotic solutions from a set of three algebraic conditions equivalent to a zero 
curvature condition, having more than one solution. It is possible, for all the models we 
describe below, to write three simpler equations, which have, as a unique solution, the 
exotic model. 

These equations define a rational variety of dimension 5 in an 8 dimensional space. They 
have special (multi)-rationality properties: for any choice of a corner in Figure [^a), the three 
corresponding variables (one vertex variable and the two adjacent bond variables) can be 
expressed rationally in terms of the five others. This special rationality property allows one 
to define rational evolutions in all four directions of the square lattice. This is a generalized 
form of the notion of quadrirationality introduced in 

Equivalently, choosing two adjacent bonds, the relations define 2 i— )■ 2 birational maps 
between the variables attached to these bonds and the two remaining ones. These maps and 
their inverses actually happen to have the same form, but they are not involutions a priori. 



The simplest case is obtained from the lattice potential KdV (HI in [T]). Computing the 
zero curvature condition on the coarser lattice we found two rational solutions [10] . We give 
here for reference the regular one, which does not fulfill condition ([T]): 
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The models 
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The exotic solution is: 

Wi2 = Wi + W2 — w, (3a) 

X2 = Y + P, (3b) 

Y, = X + P, (3c) 
_ {X -Y)[{p -r){w - W2) + {q - r){w - wi)] 



{p — r){w — W2) — {q — r){w — Wi) — {w — Wi){w — W2){X — Y) 



(3d) 



The parameters p, q, r appearing in the solutions come from the construction of the 
models: they are present in the Lax pairs we used. 

The sequence of degrees for the bond variables is 

{dn}b = 1, 4, 13, 28, 49, 76, 109, 148, 193, 244, 301 .. . (4) 
The generating function of the sequence dl]) is 

(1 - s) 

The sequence has quadratic growth and the entropy vanishes. 



Simplified form: By taking suitable linear combinations of equations we can also 
write the equations in the form 

W12 -W1-W2 + W = 0, (5a) 

X + X2-Y-Y1 = 0, (5b) 

XX2-YY,-iX-X2)^-^^ + iY~Yi)^^^ = 0. (5c) 

W — W2 W — Wi 

The equations ([5]) have the form typical for most of our results: an independent linear 
or linearizable w equation, a linear or linearizable equation for X,Y, X2,Yi, and a coupling 
equation. It is then easy to see that the evolution is rational in any direction: solve the w of 
any corner from the first equation, then its adjacent X, Y variables can be solved rationally 
from the remaining two. 



Remark: The parameters appearing in the results (here and below) follow from the Lax 
set-up |10], p,q,r corresponding to the three coordinate directions of the consistency cube. 
In the end the regular solution typically depends only on p — q, while the exotic solution 
depends on p — r and q — r. Clearly any finite r can be absorbed into p, q. The limit r — )■ 00 
is also possible, but in the present case it produces a model that is linear in X, Y . 



3.2 Hie 

This quad map was introduced in [16] . There in fact are two kinds of integrable models of 
type Hie, related by inversion. 

(x - Xi2){Xi - X2) - (P - ?)(! + ^XXi2) = (6) 

{x - xi2){xi - X2) - {p-q){l + exiX2) = 0, (7) 
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When extended to the whole lattice they must alternate thereby forming a black-white 
(checkerboard) lattice [ISl [TU] . 

We may consider two different configurations for the 2x2 basic cell 
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To each of these patterns corresponds a zero curvature condition. Each of these again 
have two rational solutions, the regular one and an exotic one. The exotic solutions are as 
follows: 



3.2.1 Hlea 

For configuration a, we have: 
Original form: 

X2 =Y + P/Q, Yi = X + P/Q, (8a) 
P ={X — Y)[{p — r){w — if2)(l + ewwi) + (g — r){w — Wi){l + eww2)\ 

+ e(p + q — 2r){j> ~ q){w — Wi){w — W2) (8b) 
Q={w-wi){w-W2){X -Y) 

— {p — r){w — 1/^2) (1 + ^wwi) + {q — r){w — Wi){\ + eww2) (8c) 
(wi + W2)(l + ew^) 

W12 + w =^—, ■ r. (8d) 

1 + e[ww\ + WW2 — W1W2) 

Entropy: The sequence of degrees are respectively: 

{dn}, = 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, . . . (9) 

for the vertex variables, and 

{dn}b = 1, 5, 17, 37, 65, 101, 145, 197, 257, 325,401, 485, 577,... (10) 

for the bond variables. The sequence ([9]) has linear growth, signaling the linearizability of the 
evolution of the vertex variables (see below). The sequence (|T0|) has the generating function 

CM = 5:d„.» = i±f^^ (11) 

1 — S 

n ^ ' 

which means quadratic growth of the degrees, i.e. integrability. 
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Simplified form: Transforming the vertex variables by the Moebius transformation w t-^ 
{w — \)/[k{w + 1)] with e = —n^ etc, takes fISdl) into 

W Wi-2, = WiW-2- (12) 

Then we can write the equations in the form 

wwi2 — WiW-2, = 0, (13a) 
X2 + X-Yi-Y = 0, (13b) 

X^X - Y,Y - (X - X,)(g - r)K^^ + (F - Y,)ip - r)n^^ 

W — W2 W — Wi 

+K^[{p-rY - {q-rf] = 0. (13c) 
In comparison with (|5]), which was linear in w, this is multiplicative. 

3.2.2 Hle/3 

For the configuration P we have the exotic solution: 
Original form: 

X2 = Y + P2/Q2, Yi = X + Pi/Qi (14a) 
={X -Y){1 + eX^)[{p -r){w- W2) + {q-r){w- Wi)], (14b) 
Qi ={X - Y){w - wi){w - W2) -{p- r){w - W2){1 + eX^) 

+ {q-r){w-Wi){l-eX{X -2Y)), (14c) 
P2 ={X - r)(l + eY^)[{p -r){w- W2) + {q - r){w - Wi)], (14d) 
Q2 ={X - Y){w - wi){w - W2) -{p- r){w - W2){1 - eY{Y - 2X)) 

+ {q-r){w-Wr){l + eY^), (14e) 
W12 =Wi + W2 — w. (14f) 



Entropy: The entropy calculation leads to the same conclusion as for the configuration a: 
the vertex evolution is linear, and independent of the bonds. The degree sequence for the 
bonds is the same as (fTU]) . 



Simplified form: Starting with (fT4|) and using the Moebius transformation X 1— j- (X 
1)/[k{X + 1)], etc. with e = — k^, we get: 



X + X2-Y -Yi + 2k— — —(X - X2) -2k— — —(Y-Yi) = 0, (15b) 

W — W2 W — Wi 



Wi2 — wi — W2 + w = 0, (15a) 

— (F-Fi) = 0, (15b) 

- Wi 

XX2-YY1 = 0. (15c) 
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Remarks: Hlea and Hie/? are deformations that in the hmit e — reduce back to HI. 
As a consequence the original forms and flUl) reduce to Q, as can be readily verified. 
However, the simplified forms f lT5]l and f[T^ are obtained with a transformation that is 
singular in e. In these equations e appears through k and the k — ?■ limit trivializes them. 

The models again depend on p — r and q — r, so r can be eliminated, except that r — )■ oo 
is a possible limit. In the Hle/3 case the r — )■ oo limit leads to a nonlinear model, which may 
be interesting on its own. 



3.3 dmKdV (HSs^o) 

If we start from the discrete modified KdV equation [H3 with 5 = in [T]) we obtain the 
following: 



Original form: 

wuw = W1W2, (16a) 

X2X = YiY, (16b) 
X2 Yi {q'^W2 — r^w){w — Wi)pX — {p'^wi — r'^w){w — W2)qY 



Y X {r'^W2 — q^w){wi — w)pY — {r'^wi — p'^w){w2 — w)qX 



(16c) 



Entropy: The sequence of degrees we get are respectively 

K}„ = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,... (17) 

for the vertex variables, and 

{dn}b = 1, 4, 13, 28, 49, 76, 109, 148, 193, 244, 301,... (18) 

for the bond variables, actually the same as for HI. The vertex evolution is linearizable, 
and the bond evolution is integrable. 



Simplified form: By taking suitable linear combinations of the original equations we can 
write the result also as 



WW12 — W1W2 = 0, (19a) 
(X + X2)p{q' + r') - (X - X2)p{q' - r'^"" ^ 



w — W2 



-{Y + Yi)q{p' + r') + {Y -Yi)q{p' -r') ■ = 0, (19b) 

W — Wi 

XX2-YY1 = 0. (19c) 



Remark: For H3 the parameter dependence is through p/r and q/r. Therefore we may 
put r = 1, but the special limits r = and r — 00 also produce interesting sub-cases. 
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3.4 dSKdV {Qls=o) 



This model was obtained from the discrete Schwarzian KdV equation [19] {Qi with 6 = 
in [1]). For the elementary square we used the Lax matrices 

L(x„x;p)=f^^^+"^^-^^) -^7^ ,V (20) 

^ ' '-^^ Y p -px + r{x - Xi)j ^ ' 

where L{xx^x\p) describes the parallel transport from x to Xi, and L{x2,x; q) from x to X2- 
When constructing the Lax pair for the bigger 2x2 lattices, we get expressions depending 
on the corner variables x, Xn, X22, 2^1122 which will be relabeled as w, Wi, W2, W12, respectively, 
while Xi, X2, X112, X122 will become bond variables X,Y,Yi, X2, respectively. 

The results also depend on the three parameters {p, q, r) in an homogeneous way. The 
"true" parameters are thus p/r and q/r. If r is not vanishing we do get a rational model, 
but that model is not integrable, as we will see. (We will not give here the explicit solution 
for r 7^ 0, as it may be extracted from the simplified form of the defining equations given at 
the end of this section.) In the r — )■ limit, we have the following defining relations: 



Original form: 

U!i2 = ^ with (21) 

£3 

A = -{X-Y){2w-X~Y)wi'^W2'^+{-2YX^-4:Yw^ + XY^ + 4:wYX + w^)wi^W2 
~2w{X - Y) {w^ - Y X) wi W2 - (-2 Y^X - 4 w'^X + AwY X + wX^ + w"^) wi w-^ 

{w - Xf -X'^iw- Yf W2^ - wY'^ {w - Xf wi + [w - Yf W2, 
B = - {2w^X -Y^X^ + AwXY^ -Aw^XY -w'^Y^)wi-w{w - Xfw2'^ 
+ (-4 w^X Y + 2Yw^- F^x^ + 4 wX'^Y - X'^w'^) W2 + w {w - Yf 
+2 (X - Y) (w^ -YX)wiW2 + {w- xf wi W2^ - {w - Yf wi'^W2 
+w'^ {X - Y) {wX + YW-2YX), 

and 

C 

X2 = — with (22) 

C = X^w -Yfw2-Y'^ {w - Xfwi + Y {w - Xfwi'^ 

- (2 Y^X - YX^ - Yw^ - 2 wYX + 2 w^X) wi W2 + {X - Y) {-X + 2w-Y) wi'^W2, 
D = - {2Y'^X -YX"^ -Yw"^ -2wYX + 2w^X)wi 

+Y {w - xf W2 + {w- Yf wi^ -{w- xf wiW2 + w{X- Y) (-2 YX + wX + wY) . 

By symmetry, Yi is obtained from X2 by the exchange (X ^ Y,wi ■v^ W2) 



Entropy: For r = the sequence of degrees for the vertex variables is 

{dn}v = 1, 10, 28, 71, 139, 248, 398, 605, 869, 1206, 1616, 2115, 2703, . . . (23) 
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fitted by the generating function 

1 + 7 s + 9 - s 
(1 + ^) (1 - s) 



u^) = : : J (24) 



For tlie bond variables we get tlie sequence of degrees 

{dn}b = 1, 8, 33, 92, 201, 376, 633, 988, 1457, 2056, 2801, 3708, 4793, . . . (25) 

fitted by 

c,(.) . '^^')';^«;^^''' (26) 

(1 - s) 

The growtli of the degrees is cubic, showing integrability. It is interesting to notice that this 
example does not have the quadratic growth so commonly observed. The vertex and bond 
variables are moreover intertwined in a non-trivial way. 

The behaviour when r = has to be contrasted with the generic r ^ behaviour, where 
the degree calculation yields for the vertex variables 

{dn}v = l, 10, 38, 149, 565, 2110, 7882, 29425, 109817, 409850, 1529590,... 

and for the bond variables 

{dn}b = 1, 8, 45, 186, 711, 2672, 9991, 37304, 139239, 519666, 1939437,... 

fitted respectively by the generating functions 



indicating non-integrability (the rate of growth of the degrees is given by the roots of — 
4s + 1, i.e. 2 + ^3). 



Simplified form: The transformation 

Xw - wi Yw-W2 ^ X2W2 - W12 ^ . YiWi - W12 

^^^3^-^^^^'^-^ X2-I -^^^ Y,-l ' 
simphfies the equations to 

X2Y-Y1X = 0, (27a) 
YiY{w-wi)+wi2-W2 = 0, (27b) 
X2X{w - W2) + Wi2 - wi = 0, (27c) 
Here is also, for reference, the (nonintegrable) form before taking the r — limit: 

X2Y-Y^X = riAX2 + BY^), 
YiY{w - wi) + W12 -W2 = 0, 
X2X(w - W2) + W12 - wi = 0, 



where 



A 
B 



{w — Wi){w — W2) qXY — ({w — W2)X — wi + W2)'^p 

pq{w - wi){wi - W2) 
{w — Wi){w — W2) pXY — {{w — Wi)Y — Wi + W2)'^q 

pq{w - W2){wi- W2) 
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4 Non-autonomous functional Yang-Baxter equations, 
three dimensional consistency 



One interesting property of the four models proposed in Sec. I3.m3.3l is that they provide 
non-autonomous solutions to the functional form of the Yang-Baxter equations [2], alias 
Yang-Baxter maps. This may be seen by examining the 3 dimensional consistency of the 
models on the configuration shown in Figure [2l The vertices introduce position dependence, 
and the maps induced on the bonds verify a modified functional Yang-Baxter equations if 
the model is 3D consistent. 



Suppose that w, wi, W2, ws, X, Y, Z are given. Then the relation on each face allows one 
to evaluate X2, X3,Yi,Y2, Zi, Z2, Wu, ^13, W23 in a unique way, X23, 113, Z12 in two different 
ways, and W123 may be computed in three different ways. Consistency means they all have 
to give the same results. We have checked that all the integrable models described here verify 
this 3 dimensional consistency. 

It should be noted that for the model presented in section 13. 4^ the consistency is not 
verified when r 7^ 0, precisely in the case where it is not integrable, but is verified for r = 
(that is to say the limit p = q = s = 00), where it is integrable, as ascertained by the 
vanishing of the algebraic entropy. 

The first four sets of simplified equations have similar structure. The vertex equation is 
independent of the bonds, and the values w might be considered as parameters. They intro- 
duce non-autonomy in the bond evolution, which turns out to be the functional equivalent 
of the quantum Yang-Baxter equation introduced in [20]. The vertex equations are either 
linear (additive case), or linearizable (multiplicative case), and can be solved explicitly. In 
the additive case w + W12 — Wi — W2 = the general solution is Wn,m = F{n) + G{m), and 
in the multiplicative case wwu — W1W2 = 0, we have Wn,m = F{n) G{m), with F and G 
arbitrary functions. The equations determining X and Y contain one linear part, and a 




W23 



Figure 2: The consistency cube 



11 



coupling equation which introduces non-hnearity. 



4.1 dpKdV and Hie, configuration a 

We can solve dpKdV ([5]) and Hie, configuration a, f|T3|l together. The vertex evolution is 
additive for ([5^) and multiplicative for flT3k) and is solved as described above. Equations 
fl5c|) and fll3cp can then be cast in the common form: 

XX2- YY, - g{m) {X - X^) + f{n) {Y-Y,)+u = 0, (28) 

where / and g are readily expressible in terms of the aforementioned F{n) and G{m) and cu 
is a constant. 

We can solve f l28|) together with the constraint X+X2 = Y+Yi for arbitrary f{n),g{m),u, 
obtaining 

X, = Y + R, Y, = X + R, with R /H'-^M'-- 



Y-X + f{n)-g{my 



Redefining the X and Y by the shifts 

X = X' + fin) + cr(n) + p(m), Y = Y' + g{m) + a{n) + p(m), 

where a, p are determined from a{n + 1) — (y{n) = 2 f{n) and p(m + 1) — p(m) = 2 g{m), we 
finally get 

xi = y' + p. y/ = x' + p, with (29) 

X — Y 

This is a non- autonomous version of the Adler map [T7| (F^ in the classification of 
4.2 Hie, configuration f5 and dmKdV {H3s^o) 



Equations f[T5|) and f[T9l) can also be solved together. In both cases the coupling equations 
fll5bl) and fll9bl) can be written as 

^{X + X2)-u{Y + Fi) - g{m) {X - X2) + f{n) {Y - Y,) = 0, (30) 

where / and g are, as above, some functions related to F{n), G{m), and p and u are constants. 
Equation fl5U|) together with the constraint X X2 — F Yi = is solved in the generic case by 

X, = YQ, Y, = XQ, (31) 

(z/ + f{n)) X - (/i + ^(m)) r 

Scaling X and F by 

X = X'a{n)p{m)/{f{n) + iy), Y = Y'a{n)p{m)/{g{m) + p), (32) 

12 



where a, p now solve 



a{n + l) = a{n)'^^jp\ p{m + 1) = p{m) ^ ^1""! , (33) 



we get the equations in the form 



with 



X'- ^' p r- ^' p p.^M^l^^MT r34i 



v'^ — fijiY^ p? — g{mY 

Equation f l5^ is a non-autonomous version of F/// of 

It should be noticed that the non-autonomous nature of the above equations parallels 
the results of 



(35a) 



4.3 dSKdV (gij^o) 

In the case of equation (127|) the w evolution is not linearizable, as indicated by the algebraic 
entropy analysis. If we consider (1271) as three equations for two variables X, Y then we can 
derive from their compatibility an equation for w, which can be written as 

T(n,m + 2) _ H{n + 2,m) 
T{n,m) H{n,m) ' 

where T, H are 3-point Schwarzian-like derivatives of n, m, respectively: 

. w(n + 2,m) - w(n + l,m) , w(n,m + 2) - w(n,m + 1) , 

T(n,m):= -— — , H{n,m) := -— — . (35b) 

w[n + l,m) — w[n, m) w[n, m + 1) — w[n, m) 

Equation ([35]) connects the points of a 3 x 3 sublattice, except the center point (in the above 
w(n -|- 1, m + 1)). The algebraic entropy analysis indicates the rare cubic growth. Equations 
defined on a 3 x 3 sublattice are typical for Boussinesq-type lattice equations [22], but they 
usually have quadratic growth. 

The Schwarzian KdV equation is given by 

, [w(n, m) — w(?T. -|- 1, m)][w(n, m + 1) — w(n + 1, m + 1)] p 

[w[n, m) — w[n,m + 1)\ [w[n + 1, mj — w{n + 1, m + IjJ q 



From fl35l) one can now derive 

S{n + 1, m + l)S{n, m) = S{n + 1, m)S{n, m + 1), (37) 
which can be solved and we find that w solves the non-autonomous SKdV equation 
[w(n, m) — win + 1, m)] [w{n, m+1) — w{n + 1, m + 1)] f{n) 



[w{n, m) — w{n, m + 1)] [w{n + 1, m) — 'w{n -|- 1, m + 1)] gifn) 



(38) 



Thus when the model (1271) is considered on larger cells made of four adjacent squares, we 
find an independent non-trivial evolution of the vertex variables (135|) . driving the evolution 
of the bond variables according to (I27b .c). 
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5 Conclusion and perspectives 



We have presented lattice models with values given at both the vertices and the bonds of a 
2D square lattice. They were constructed using the ambiguity of the zero curvature condition 
on a coarse grained lattice. The integrability of these models is guaranteed by the vanishing 
of the algebraic entropy and by the Consistency- Around-the-Cube property. 

Four out of five models provide non autonomous generalizations of known Yang-Baxter 
maps, the fifth case (1271) (related to dSKdV) being seemingly different. 

They all share some remarkable algebraic properties, in particular their multi-rationality, 
which might be the frame for further examples and an eventual classification. 

Many more aspects will have to be examined (especially for (1271) ). like the existence 
of Backlund transforms, proper Lax pairs, symmetries, reductions (periodic reductions as 
well as similarity reductions) which should produce interesting integrable maps, including 
discrete Painleve equations. Their continuous limits should also be considered. All these go 
beyond the scope of this paper, and require further studies. 
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